General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




Electro-Optics Systems Laboratory 

Department of Electrical and Computer Engineering 



(HASA-CB-176026) PHASE ICCKEE LOOP M65-3C178 

SmCHa0Nl21II0H FOB DIBECI DE1EC1IOM OPIICAI 
PPH C OH HU IIC ATX C N SYSTEHS (Illinois Unit., 

Urban a-Cha«paign« ) 53 p HC A04/HF A01 Unclas 

CSC I 17B G3/32 21785 




Phase-Locked Loop Synchronization for Direct 
Detection Optical PPM Communication Systems 



University of Illinois at Urbana-Champaign 





UILU-ENG-85-2549 


PHASE LOCKED LOOP SYNCHRONIZATION FOR 
DIRECT DETECTION OPTICAL PPM COMMUNICATION SYSTEMS 


by 

C. C. Chen 
C. S. Gardner 


EOSL Publication No. 85-003 


Technical Report 
May 1985 


Supported by 

Contract No. NASA NSG-5049 


NATIONAL AERONAUTICS & SPACE ADMINISTRATION 
Goddard Space Flight Center 
Greenbelt, Maryland 2G771 


ELECTRO-OPTIC SYSTEMS LABORATORY 
DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING 
COLLEGE OF ENGINEERING 
UNIVERSITY OF ILLINOIS 
URBANA, ILLINOIS 61801 



ABSTRACT 


Receiver timing synchronization of an optical PPM communication system 
can be achieved using a phase-locked loop (PLL) if the photodetector output is 
properly processed. The synchronization performance is shown to improve with 
increasing signal power and decreasing loop bandwidth. Bit error rate (BER) of 
the PLL synchronized PPM system is analyzed and compared to that for the 
perfectly synchronized system. It is shown that the increase in signal power 
needed to compensate for the imperfect synchronization is small (less than 0.1 
dB) for loop bandwidths less than 0.1 Z of the slot frequency. 
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1 . INTRODUCTION 


M-ary Pulse Position Modulation (M-PPM) has been shown to be an effective 
modulation technique for direct detection optical communications. In M-ary PPM 
systems, each word frame is divided into M time slots and the data are encoded by 
transmitting a single laser pulse during one cf the M time slots. The perfor- 
mance of optical PPM systems has been well documented [1,2,3]. However, most 
studies were based on the assumption of perfect timing synchronization between 
the receiver and the transmitter. For a practical communication system, this 
condition is not always satisfied. 

In the presence of timing error, some of the signal photons that are trans- 
mitted in the signal time slot will be counted in the adjacent slots. This 
effectively reduces the signal-to-noise ratio and, consequently, increases the 
probability of decoding error [3]. Accurate timing synchronization is needed 
to minimize the decoding error at the receiver. In general, timing synchroni- 
zation can be achieved by either transmitting a separate timing signal along 
with the data, or by acquiring synchronization directly from the received signal. 
For applications where the transmitter power is limited, the latter method is 
preferred. Phase-locked loops (PLLs) provide an easily implemented method for 
recovering the transmitter timing. The applications of phase-locked loops in 
radio frequency systems have been studied extensively [4,5]. The use of phase- 
locked loops for timing synchronization in optical communication systems have 
also been studied in recent years [6,7,8]. Gagliardi and Haney [6] and Snyder 
and *. arrester [7] analyzed the probability density of tracking error for a PLL 
under shot noise input. Mengali and Pezzani [8] studied the phase error variance 
of a PLL driven by photodetector current in an optical pulse amplitude modulation 
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(PAM) system. These studies have shown that that PLLs can be used to track the 
transmitter timing provided that the transmitted signal contains a frequency 
component at the desired lock-in frequency. In related work, Georghiades [9] 
and Georghiades and Snyder [10] studied the problem of obtaining frame and word 
synchronizations of a PPM system using a coding technique, with the assumption 
that slot synchronization has been achieved. 

In this paper, It is shown that for a PPM system which transmits square 
pulses occupying the entire signal slot, pha3e-locked loops cannot lock onto the 
photodetector output directly because the detected signal does uot contain a 
spectral component at the time slot frequency. However, with proper preproc- 
essing of the photodetector output, namely, squaring the detected signal, PLLs 
can be used to lock onto the transmitter slot frequency. The performance of the 
PLL using the preprocessed PPM signal is then analyzed using a perturbation 
method [8], and an expression for phase error variance is derived. Finally, the 
error performances of phase-locked loop synchronized PPM systems are evaluated and 
compared with perfectly synchronized systems. 
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2. SIGNAL AND NOISE AT THE PHOTODETECTOR OUTPUT 

The output current <f the photodetector can be modeled as the sum of a signal 
shot noise and a Causslan thermal noise. Furthermore, the signal shot noise can 
be modeled as a filtered point process in which the output is the superposition 
of the detector response to each detected photon [11]. The detector output can 
be written as 

«t) - yo + i th (t) - 1 g . h<t - Tj ) + i th (t) (i) 

'V 

where i (t) and i . (t) are the signal shot noise and the thermal noise current, 
s th 

respectively, h(t) is the combined pulse response of the detector and subsequent 
amplifier-filters, G is the photodetector gain (assumed to be constant), and 
is the arrival time of the j C ^ photon. For simplicity, it is assumed that the 
internal gain of the photodetector is large so that the effects of the thermal 
noise can be neglected. The principal source of noise at the detector output is 
therefore the signal shot noise, which is due to the stochastic nature of the 
photon counting process. 

The photocount statistics at the photodetector, conditioned on the received 
optical power, can be shown to be Poisson distributed with count rate X(t), 
which is related to the received optical power [2]. 

X(t) = qP (t)/hv . (2) 

o o 

n is the quantum efficiency of the detector, P Q (t) is the total received power 

at the detector surface and hv is the photon energy. The total optical power 

o 

incident on the detector surface is the sum of the received signal pover and 
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the power of the background radiation. The photocount rate A(t) is therefore 
the sum of the background rate and the count rate due to the signal. For an 
M-ary PPM system, the photocounc rate can be modeled as 

X(t) - \ * l X, p(t - ki w - y,) (3) 

k 

where A is the photocount rate due to the signal, A. is the constant photocour.L 
s d 

rate due to the uniform background radiation, T is the word trame period, T is 

w s 

the time slot width, which equals T^/M, is codeword which takes on 

the integer values (0,1 M-l) , and p(t) is a unit square pulse of width T g . 

Since the transmitted data {C^} are random, A(t) is a stochastic process. 
The output of the photodetector can therefore be regarded as a sample function 
of a doubly stochastic filtered Poisson process [12]. The statistics of this 
process can be evaluated by first taking the expectation conditioned on the 
received photocount rate A(t), using the fact that the conditioned process is a 
filtered Poisson process. This conditioning can then be removed by taking the 
expectation with respect to the transmitted data. The expectations of the 
filtered Poisson process can be obtained by differentiating its characteristic 
function [12], given by 

*(w) - E[e" Jwi s (t) ] 

- exp[ / AC C ) ( e ^-l)d£] . (4) 

The first and second moments of the detector output, conditioned on the 
detector photocount rate, are given by 
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E[i(t)|x] 3 GX(f *h(t) 

(5) 

E[i I 2 (t)|x] 3 G 2 [X(t)*h 2 (t) + (X(t)*h(t)) 2 ] . 

(6) 


The unconditioned moments of the detector output can be obtained by taking the 
expectations of Equations (5) and (6) with respect to the data {C^}. If {C^| 
are independent and uniformly distributed, the resulting unconditioned moments 
of the detector output are given by 

K m 

E[i(t)] =■ G(X b + ^ l p(t - lT s ))*h(t) (7) 

X 00 X 00 

E[i 2 (C) ] - G 2 (U b + jp l P(t - IT )}*h 2 (t) + 1(\ + jp l P(t - lT s ))*h(t)] 2 

00 l=-co 

I O ® . . co M-l 

+£x 2 I [h(t)*p(t - iT g )] Z ” ^2 l l h(t)*p(t - kT w - vT s ) 

£*-oo M k=»-ao y=0 

V=0 

X h(t)*p(t - kT^ - y Tg ) ] (8) 

where p(t)*h(t) denotes the convolution of p(t) and h(t). Notice that by taking 
the expectation with respect to the transmitted data, the expected moments of 
the detector output now contain periodic terms with period t • Phase-locked 
loops can then be used to lock onto these signals and generate timing references. 
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3. SIGNAL PREPROCESSING 


In an optical M-ary PPM signaling system that transmits equally likely 
codewords, the optimal receiver [1,2] compares the received photocount in each 
of tne M time slots and chooses the time slot with the largest photocount. This 
decoding operation cannot be achieved without timing synchronization between the 
receiver and the transmitter. For this reason, i. is important for the receiver 
to maintain a local timing reference which identifies hue beginning of each 
transmitted time slot. With proper preprocessing, phase-locked loops can be 
used to generate this local timing reference from the photodetector output. 

A block diagram of a typical phase-locked loop [4,5] is shown in Figure 1. 
It consists of a voltage controlled oscillator (VCO), a low pass filter and a 
phase detector, which generates an output voltage as a function of the phase 
difference between the VCO output and the incoming signal. In order to rchieve 
effective phase synchronization, the input signal to the PLL must contain a 
strong periodic component at the lock-in frequency. 


From Equation (7), it is seen that if the transmitted pulse p(t) is a square 


pulse of width T g , the expected output of the photodetector is a constant 

X 00 

E[i(t) ] - (X b + jj— )*h(t) - (X b + m“)*/ h (f)dT . 


(9) 


Therefore, the detector output does not contain a periodic comDonent at the 
slot frequency. This fact is also seen by examining the signal power spectrum. 
Based on the model of the photodetector output given in Equation (l), the power 
spectrum of the detector output for an M-ary PPM receiver can be shown to be 
(Appendix A) : 
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S(w) 



|HU)| 2 .[(X b + yp) + 2ir5(u)(X b 2 
lp(a))| 2 (l - |R( w) | 2 ) + ~ k b 


2X sV 

M ; 


| pv uj) |‘ 


• \ 5( a) 

k — ® 



( 10 ) 


where 


R(u) 


E[exp(-jC^aT g ) ] 


M-l 

ff I 

t-0 


“JWllTg 

e 


(11) 


and P(w) and H(uj) denote the Fourier transforms of the pulse shape and the 

impulse response of the detector-filter, respectively. If the transmitted pulse 

shape is a square pulse with width T , its Fourier transform P( w) will be zero 

s 

at all integer multiples of the slot frequency, 2 tt/T . In this case, it Is seen 

s 

from Equation (10) that no discrete frequency components at the slot frequency 
or its harmonics exist in the detector output spectrum. Consequently, the phase- 
locked loop cannot track the output of the photodetector directly. Preprocessing 
of the detector output is necessary to generate a frequency component at the 
slot frequency. 


Only nonlinear processing needs to be considered, because any linear filter- 
ing will simply amount to multiplying the power spectrum in Equation (10) by the 
magnitude square of the filter transfer function, and will not affect the absence 
of the spectral component at the slot frequency. One approach is to filter the 
signal and then square the filter output as depic ’ in Figure 2. The output 
of this preprocessing circuit can be written as 


i 2 (t) - i 2 (t) 


lT j 


G*h( t 

} 



2 


( 12 ) 
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Figure 2. The structure of the receiver timing subsection. The output of 
the photodetector is firr»t filtered by a preprocessing filter 
before being squared. Tb** output of the squarer la then fed 
into the PLL to genera-. r \ v ocal timing reference. 







where h(c) now denotes the combined Impulse response of the detector and the 
preprocessing filter. The expected output of this preprocessing circuit Is 
given by Equation (8). By examining Equation (8), It Is seen that the third term 

at the right-handed side is periodic with period T and, unlike the photodetector 

s 

output, does not liave a vanishing Fourier transform at the slot frequency. The 
output of the preprocessing circuit can therefore be regarded as the sum of a 
periodic signal s(t) and an additive noise term n(t), related by 

i 2 (t) - s(t) + n(c) (13) 

where from Equation (8) 

•' X a 2 " 2 
s( t) * G* -g- l [p(t)*h(t - *T g )r 

t =-«0 

n(t) - i 2 (t) - s ( t ) . (14) 

The power spectrum of the preprocessed signal is quite complicated. However, 
the expression for the power spectrum can be simplified considerably if it is 
assumed that H(a)) blocks the dc component of the signal, and that both P(w) and 
H(u>), the Fourier transforms of p(t) and h(t), are slowly varying functions of 
frequency compared to R(u>), defined in Equation (11). The first assumption is 
made because the dc component of the detector output contains no timing informa- 
tion and will only contribute noise to the squared signal. Therefore, it is 
desirable to remove the dc component from the photodeteccor output. The second 
assumption holds for higher-order PPMs. Because, in general, h(t) and p(t) are 
pulses of width comparable to T , their Fourier transforms will have supports on 

the order of hi j 2ir/T . The support for Run), on the other hand, is on the 
s s 



order of w g /M, where M is the PPM order. Therefore, for high-order PPM, the 
assumption that P(w) and H(w) are slowly varying compared to R(w) usually holds. 
For high signal counts, the power spectrum of the preprocessed signal 1^(0 near 
the slot frequency can be approximated under these assumptions by (Appendix B) 


S (w) - S ( u»> + S (u) 

Z s n 

whera 

x 4 

A / I n 

S (dJ - |P(to)H(aj)*r(a))H((o) | *6(to - to g ) 

3 

S < u ) - 4 _s_G 4 . 1 f ?(a» , )P(u")P(a)' - w" )H( ,V )H( a> - tt f )H?u") 

n T v2 

w (2ir) 


* 

x H( to — uDdto'duf . 


05) 


(16) 


07) 


The power spectrum of the preprocessed photodetector output car be separated 
#a*o two terms. S g (w) is due to the sinusoidal signal component s(t) in the 
preprocessed signal, and S n (u)) is the power spectrum of the additive noise, 
n(t), which is the random component of the preprocessor output. From Equation 
(16), s(t) can be written a3 

8( t) - A sin(iu t + 0 ) 
s s 

n 2 
2X G 

A ■ — ® I •=— f P ( ti3 )H(co* )P(to “ to')H(aj - <*>* )du>* | (18) 

i * Zir J s s 

w 


where u g = 2ir/T g is the frequency of the sinusoid, which is equal to the slot 
frequency, and A and 0^ are the amplitude and phase, respectively. 



4. PLL ANALYSIS 


The presence of the sinusoidal corap lent with slot frequency at the preproc- 
essor output indicates that PLL can now be used to track the transmitting timing. 
The PPL depicted in Figure 1 has the equivalent baseband model shown in Figure 3 
[4,5]. In reducing the PLL to its baseband model, it is assumed that the input 
consists of a sinusoidal signal s(t) with frequency w g and amplitude A, and an 
additive noise term n(t), and that the V C0 has gain constant and free-running 
frequency Wq, which differs from w g by an amount ft, the frequency detuning. 
From the baseband model, it is straightforward to show that the phase error $ 
between the VCO output and the incoming sinusoid satisfies 

"ft - K^f(t)*[-j sin$ + n(t) cos(w g t + 9 g “ $)] (19) 


where f(t) is the impulse response of the loop filter, and 9 is the initial 

s 

phase of the incoming signal. In general, Equation (19) is a stochastic integro- 
differential equation that is difficult to solve. However, if the noise term 
n(t) is small, we can make a perturbation expansion of Equation (19) [8]. The 
result, after expansion, is a series of equations that can be solved recursively. 
The first two terms of the expansion satisfy the equations 


d $0 A 

d — " " • 2 K v f(O * sln4, 0 
d $ 1 A 

— = -K v f(t)*(- 2 * 1 cos 4 > 0 + n(t)cos(u) s t + 9 g - <fr Q ) ) • 


( 20 ) 

( 21 ) 


The phase error at the VCO output can then be approximated by 
$(t) = $ Q (t) + ^(t) . 


( 22 ) 



n(t)Cos(u> s t +9-<f>) 




= & - KF(p) { A Sin <£ + n(t)Cos(w s f +9 -<£)} 


Figure 3. Baseband model of the phase-locked loop. 





The solution of Equation (20) gives the phase error as a function of time 
when no noise appears at the input. In the steady state, <J>q tends to 

<fr n (») » sin *(- — ) (23) 

U y KF(0) 

where F(s) is the transfer function of the loop filter, which is related to the 
impulse response f(t) by a Laplace transform. When noise is present at the 
input to the PLL, higher-order solutions to the phase error must be included to 
account for the noise effects. The first-order equation, Equation (21), contains 
a stochastic driving term n(t) on the right-hand side. Consequently, its solu- 
tion will also be stochastic. In the steady state, the solution to Equation (21) 
is given by 

♦ j(t) a g(t)*[n(t)cos(uj g t + 9 s - 4 >q) 1 ( 2/ 0 

where g(t) is the inverse Laplace transform of the loop transfer function G(s) 
defined as 

G(s) = j KF(s)/(s + KF(s)cos$ 0 ) . (25) 

Equation (24) states that the first-order solution of the phase error is the 
filtered output of the stochastic driving term 

v(t) - n( t)cos( uj t + 0 - $ ) . (26) 

s s U 

The autocorrelation function of the phase error R (t,r) can therefore be related 

4 > 

to the autocorrelation of the process v(t) by 
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from which the variance of phase error can be calculated. In general, $(t) is 

not stationary. In fact, it can be shown that for PPM transmission, n(t), and 

consequently v(t), are cyclostationary with period T • It follows that the vari- 

2 

ance of phase error is also periodic with period T^. 

However, when the loop bandwidth is small, it is found that $ is actually a 
wide sense stationary (WSS) process [13]. More specifically, $ is WSS if the 
loop transfer function satisfies 

|G(u») I - 0 for OJ ; Y" = ^ (28) 

w 

In this case, the variance of $ can be written as 

a ^ * 4 — f | G( a>) | ^ • S~(u>)duj (29) 

$ ** pf v 

where S^(ai) is the power spectral density of the WSS process v(t) derived from 
v(t) by randomizing its initial phase, i.e., letting v(t) = v(t + x) with x, the 

initial phase, to be uniformly distributed over (0,1^). In practice, for the 

PLL to have a sufficiently small phase error, the loop bandwidth should be much 
smaller than the time slot frequency. Consequently, Equation (28) is almost 
always satisfied. 


The power spectrum S~(<d) can be evaluated in terms of the noise power spec- 


trum, which yields 

S~(w) ■ 7 - (S (u» + at ) + S (uj - a) ) 1 
'* 4 v n s n s J 


+ £im -4= / e / R (t,t + t)cos(u) (2t + t) + 26 - ^.Idtdi (30) 

T 2T J in v s s 0' 

I -►« —» — T 
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where S (ui) and R (t,x) are the time-averaged power spectrum (Appendix A) and 
i. n 

the autocorrelation function of the noise, respectively. Fcr the PPM signaling 
scheme, the last term In Equation (30) is negligible compared to the first two 
terms. The variance of phase error in Equation (29) can then be written as 



^2 hi I l G(t,,) | 2 ( S n (a) + u s ) + S n (u ~ <0 s ) ) dw * 


(31) 


Equation (31) can be further simplified if the bandwidth of the loop transfer 

function is sufficiently small so that S n (w) is approximately constant over the 

support of G(oj). Using the fact that S (-to) ■* S (w), Equation (31) can be 

n n 

approximated by 


b l ( ,1 

A 2 , " 8 p 


(32) 


where p = A n/fS (u: )B, ) can be interpreted as the signal-to-noise ratio of the 
v n s L' 

synchronization signal within the effective loop bandwidth B (radians/second), 

la 

defined as 


O 

b l t I |G(w)| dm. 


(33) 


By substituting the expression for noise power spectrum from Equation (17) 
and the expression for signal amplitude from Equation (18) into Equation (32), 
the following expression for the variance of phase error is obtained: 


2 

% * (T^ y 

r s s 


( 34 ) 


[ 


2T jj P((j)’ )P( w”)ptw' ~ u)")H(oj - a)' )Hta)")H( id - uDdaj'duj" 


|P(aj)H(a))*P(aj)H(aj) 


a) ■ u 



where K is the signal count per word, M is the order of the PPM, B is the 
equivalent loop bandwidth, oj g is the slot frequency, and y is a dimensionless 
parameter, which depends only on the pulse shape and the preprocessing filter 
transfer function. The values of y for some choices of preprocessing filter 
transfer functions are listed in Table 1. For a given pulse chape, the transfer 
functiou of the preprocessing filter can be chosen to minimize the value of y. 
One choice is to model H(co) as an ideal low pass differentiator with bandwidth 
Bq. Figure 4 is a plot of the value of y versus the bandwidth of the differenti- 
ator for the case where p(t) is a square pulse of width T g . It is shown that 

the value of y is minimized (y * 3.2) for E = 1.3w . Also shown in Figure 4 is 

u s 

the value of y evaluated using the preprocessing filter which consists of a low 
pass Gaussian filter with rms bandwidth Bq, followed by an ideal differentiator. 
It is seen that the minimum values of y obtained using both preprocessing filters 
are similar. 

Equation (34) shows that variance of the phase error is inversely propor- 
tional to the signal power (i.e., signal photocount K g ) and is proportional to 
the number of time slots M and the loop bandwidth. It should be noted that the 
loop bandwidth B , defined in Equation (33), actually increases with increasing 

L 

signal amplitude. In order to accommodate a wide dynamic range of input signals, 
it is the usual practice to precede the PLL by either an automatic gain control 
(AGC) circuit or a limiter. The effect of the AGC or limiter is to dynamically 
scale the input signal so that its amplitude remains essentially constant and, 
therefore, the loop bandwidth also remains constant. By using an AGC or limiter, 
it is seen that the performance of the PLL will improve inversely with the 


signal strength. 
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0.40 0.80 1.20 1.60 2.00 


DIFFERENTIATOR BANDWIDTH (B 0 /<j s ) 


Figure 4. The value of y as a function of the differentiator bandwidth for 
PPM systems with square transmitted pulses of width T g . H( m) is 
modeled as an ideal low pass differentiator, where Bq is the- 
cutoff frequency, and as a Gaussian low pass filter, where Bg is 
the RMS bandwidth, followed by a differentiator. 
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5. COMPUTER SIMULATION 

As was shown in the previous section, the timing reference generated by the 

2 

PLL synchronizing circuit is in general a random variable with variance o , 

<P 

which is inversely proportional to the signal photocount, and is proportional to 
the order of PPM and the loop bandwidth. 

In order to verify the expression of the variance of the phase error, a 
computer simulation of the phase-locked loop was developed. First, the output 
of the photodetector was generated as a Poisson arrival process with the count 
rate given by Equation (3). This photodetector output was then filtered and 
squared to simulate the effect of the preprocessing circuit. The output of the 
preprocessing circuit was then fed to the phase-locked loop simulator written 
in Advanced Continuous Simulation Language (ACSL) running on a CYBER- 175. The 
simulator, modeled after the PLL in Figure 1, had the characteristics listed in 
Table 2. 

Several simulations were carried out for different signal levels, each for 
a time period equal to 10,000 time slots, with the effective loop bandwidth, 
B^, held constant for each simulation by adjusting the photodetector gain. The 
simulator calculated and recorded the phase error at the voltage controlled 
oscillator at the beginning of each time slot. The phase error variance of the 
VCO was then calculated from the simulated data. The results of the simulation 
are shown in Figures 5 through 7. 

Figure 5a shows the power spectrum of the photodetector output before pre- 
processing. Notice the absence of the spectral peak at the slot frequency. 



Table 2 


Characteristics of the PLL Simulator 


Normalized Signal 
Amplitude 

A -i 

Loop Filter 

F(s) 20s + 1 

VCO Gain Constant 

K v = 0.33, 0.04 

Loop Bandwidth 

®L a 1 

*» AK 

<i) 2 v 

s 

Photodetector Gain 

o ttT , 

2 w 

x 2 |p( io)H( o,)*p( oj)H( oi) | 

S (j) ” Ug 

Preprocessing Filter 

h(t) = r sin t 0 < t < T s 

1 0 s elsewhere 

Background Signal 
Level 

o 

n 

PQ 













Figure 5fc shows the spectrum of the preprocessed signal in which the spectral 
component at the slot frequency is clearly seen. Figures 6 and 7 contain olots 
of the variance of the phase error versus the signal count rate and the PPM 
order, respectively. The results show that the variance of phase error is pro- 
portional to the order of the PPM and is inversely proportional to the signal 
count rate. These results are in general agreement with Equation (34). 



POWER SPECTRUM g POWER SPECTRUM 



0.00 1.00 2.00 3.00 4.00 

FREQUENCY (w/w s > 


5a. Power spectrum of the simulator photodetactor output for a 4-ary 
PPM system whici. transmits square pu?ses occupying the. entire 
signal slot. Nc spectral peak is present at the slot frequency. 



0.00 1.00 2.00 3.00 4.00 

FREQUENCY (cu/c j s ) 


Figure 5b. Power spectrum of the preprocessed simulated photodetector out- 
put. A spectral peak at the slot frequency is clearly evident. 
The preprocessing filter impulse response h(t) used in this 
simulation was 

u/f.'v m t sin(u)gt) 0 < t < T s , 

' ' i 0 elsewhere. 




PHASE ERROR VARIANCE (rod 



Figure 6. Computer simulated phase error variance versus the signal count 
rate for a 4-ary PPM receiver. The loop bandwidth and the order 
of the PPM are held constant during the simulation, and the VCO 
gain constant Ky is set to 0.33. The solid curve is the theore- 
tical variance predicted by Equation (34). 




PHASE ERROR VARIANCE (rad 11 ) 
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Figure 7. Computer simulated phase error variance versus the order of PPM 
for an M-ary PPM receiver. The loop bandwidth is held constant 
during the simulation by adjusting the gain of the preprocessing 
filter, and Ky is set to 0.04. The resulting phase error vari- 
ance is proportional to the order of the PPM. The data points 
are the results of the simulation, and the solid line is the 
theoretical variance predicted by Equation (34). 
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6. PERFORMANCE OF PLL SYNCHRONIZED M-ARY PPM SYSTEM 


The probability of word error (PWE) of * perfectly synchronized, shot noise 
Halted M-ary PPM receiver using a photon counting detector has been shown to be 


[ 1 , 2 ] 

PWE(K K^.H) - 1 -^’ (K S + *b> - (j £j<K + • e- (K s + K b ) 

k-1 



£-0 £! 


(35) 


where K g - A g T g a nd * X^T g are the expected photocounts per slot due to the 
signal and background radiation respectively. The PWE given in Equation (35) 
is difficult to evaluate because of the large number of summations. For PPM 
systems with a large signal-to-nolse ratio, the union bound can be used to give 
a good approximation of the PWE, 


PWE(K g ,K b ,M) < (M-l).PWE(K g ,K b ,2) . (36) 

PWE(K g ,K^,2) is the PWE of the binary PPM system, which is given by [14] 

PWE(K g ,K b ,2) £ y (1 + Q(/2K^, /2ip - Q(/2 T*. /2K^) ) (37) 

where Q(a,B) is Marcum's Q function, defined by 

Q(a,B) - /Y (a2 + %1)n 1 0 ( ® ),xdx ' 

3 


( 38 ) 
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The union bound of the PWE for a 4-ary PPM system is plotted along with the 
exact error rate evaluated by Equation (35) in Figure 8. The result shows an 
excellent agreement for large signal levels. 

In the presence of receiver timing error, signal photons that are transmitted 
during the signal interval may he counted in the adjacent background slots, 
thereby increasing the probability of decoding error. For a fixed timing offset 
AT at the receiver, the expression for the PWE is very complicated and difficult 
to evaulate. However, if the amount of timing error is small, the union bound 
can again be used to aporoximate the PWE. Since the effects of the fixed timing 
error are to decrease th»* expected photocount in the signal slot to increase 
the expected photocount in one of the two adjacent slots, we can wr^te 


PWE(K ,K ,M; e ) * PWE(K ' ,K ',2) + (M-2)PWE(K ' ,K ,2) (39) 

s b s b s b 

where we have expressed the timing offset AT in terms of che normalized timing 

error e = AT/T g , and K g * » ((1 - e)X g + ^) T S an< * * (eX g + X^)T g are expected 

counts in the signal and affected adjacent slots, respectively. The first term 

on the right-hand side of Equation (39) is the error rate of the binary PPM 

system consisting of the signal time slot and the "contaminated” background slot, 

and the second term is the error rate of the (M-2) binary PPM systems composed 

of the signal slot and one of the remaining background slots. Figure 9 is a 

plot of the PWE of a 4-ary PPM system versus the fixed riming error, e. The 

degradation of PWE with increasing e is clearly seen. 


Equation (39) gives an upper bound for PWE in the presence of a fixed timing 
error. For receivers employing dynamic phase synchronization circuits such as 
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SIGNAL COUNT / WORD (K s ) 


Figure 8. Comparison of the exact error rate and the union bound estimate 
for a 4-ary PFM system. Broken curves are the union bound given 
by Equation (36) and solid curves are the exact error rate given 
by Equation (35). 
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a PLL, the receiver timing e.ror would In general be random. The exact proba- 
bility distribution of the timing error Is complicated to derive. From Equation 
(24), it Is seen that the exact distribution of the timing error depends on the 
detailed statistical properties of the noise. However, when the loop bandwidth 
is small compared to the word frequency, the phase error can be regarded as a 
weighed average of the noise v(t), given In Equation (26), over many periods of 
the transmitted word. Because v(t) is essentially uncorrelated over different 
transmitted words, it is seen that the probability distribution of the timing 
error will be approximately Gaussian. A similar problem was studied by Gagliardi 
and Haney [6], in which they show that the phase distribution for a ’’LL driven 
by a shot noise limited synchronization channel was given by 

v*> - e 2?y c :r' ■!♦!<« (4o) 

where 1^ is the modified Bessel function of order zero, and a is the average 

number of sync-signal counts occurring in the time period 1/2B L » For a well- 

synchronized system, <} is small, and P can again be approximated by a Gaussian. 

♦ 

For a PIL synchronized system, the timing error renaius essentially constant 
over a given word because the loop bandwidth B^ is usually much smaller than the 
word frequency- The unconditioned error probability for such a system is there- 
fore the expectation of Equation (39) with respect to the distribution of e, which 
is approximately Gaussian, 

, 2,-2 

- = - / PVE(K ,tt ,M,e)e' E /2<I e de 

/2wct s 15 

e 


PHE< W M) 


(41) 
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Equation (41) is difficult to evaluate in a closed form. In general, numerical 

integration is required. Figure 10 is a plot of the PWE versus the expected 

signal count, K , for a 4-ary PPM system in which the receiver timing error is 
8 

assumed to be Gaussian distributed. The result shows a severe degradation of 
receiver performance with increasing timing error variance. 


The timing error variance o for a PLL synchronized M-ary PPM system can be 

e 


related to the phase error variance o derived in the previous section by 

♦ 



o 2 /( 2 tt) 2 
9 


< 2») 2 


YMB 

L 


K m 
s s 


<at 2 > 



s 


(42) 


Combining (41) and (42), the PWE of a PLL synchronized PPM system can now be 
expressed as a function of K g , K^, M, B^» and y* The effect of these parameters 
on the performance of the PLL synchronized PPM system can then be studied and 
compared to the perfectly synchronized system. 


When comparing the performance of different communication systems, it is 
usually desirable to express the error probability in terms of the probability 
of bit error (PBE). For M-ary PPM, the PBE is related to the PWE by 

PB ‘ ? * “ { * PWE * (43) 

Numerical evaluation of PBE has been carried out for various signaling condi- 
tions and receiver loop bandwidths. The results are shown in Figures 11 through 
14. Figure 11 is a plot o' the PBE versus signal count for various loop band- 

_3 

widths. It is seen that for small loop bandwidths (B / u> <10 ), the performance 

L S 

of the PLL synchronized system is almost identical to that of the perfectly 
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SIGNAL COUNT /WORD (K s ) 

Figure .0. Probability of word error (PWE) of a 4-ary PPM receiver as a 
function of K^, the expected signal photocount per word, for 
different values of the receiver timing error variance. The 
timing error is assumed to be Gaussian distributed. Background 
count rate - 5. 
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synchronized system. For large loop bandwldths (B / <*> >10 ), however, higher 

L s 

signal levels are needed to compensate for the effects of synchronization errors. 

Figure 12 is a plot of the PBE versus signal count for various background noise 

levels and loop bandwldths. Again, It Is seen that when the loop bandwidth Is 

-3 

sufficiently small (B /w <10 ), only small Increases In signal power are needed 

to achieve the same performance as the perfectly synchronized system. The 

Increases In signal power needed to compensate for the imperfect synchronization 

can be described in terms of a system loss factor. Figure 13 is a plot of the 

-9 

loss factor versus the loop bandwidth at a fixed PBE of 10 for a 4-ary and 

a 8-ary PPM system. It is seen that for small loop bandwldths (B / co <10 ^) , 

Li S 

the imperfect synchronization accounts for only 0.1 decibel loss in system 
performance, while at higher loop bandwldths ( B / uj = 10 *), the loss can be 

Lt S 

significant. 

The effect of increasing PPM order on the performance of the PPM system was 

also studied. The PBE of a perfectly synchronized PPM system decreases with 

increasing PPM order. For PPL synchronized PPM systems, however, the variat^.i 

of PBE with PPM order is more complicated, because the timing error is also a 

function of the PPM order. Figure 14 is a plot of the PBE versus PPM order for 

various loop bandwldths and background count rates with the number of signal 

photons per bit and the equivalent bit period kept constant. It shows that for 

-4 

small loop bandwidths (B / uj * 10 ), the performance of the PLL synchronized 

PPM system is almost indistinguishable from the perfectly synchronized system at 
low PPM orders (M<10). At higher PPM orders, however, because -he phase error 
variance increases with PPM order, the degradation in performance will become 
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Figure 11. Probability of bit error (PBE) of a PLL synchronized 4-ary PPM 
system versus signal count/word fo. various loop bandwidths. 
10 






LOSS FACTOR (dB) 


Figure 13. 



dB loss of transmitter energy versus the loop ban..«/J dth for i 
PLL synchronized PPM systeir at PBE = 10“9 • Kjj = 1 , y - 10. 
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more apparent. In addition, since the synchronization error increases with PPM 
order, It is expected that the PBE will become an increasing function of PPM 
order at higher order PPM. Consequently, an optimal order of PPM exists that 
will minimize the PBE. Figure 14 shows that for a PPM system operating with a 
received signal level of two photons per bit and a background count rate of one 
photon ber bit interval, the PBE is minimized at M = 100. 
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7. CONCLUSION 

Phase-locked loops can be ust to synchronize slot timing for an optical PPM 
system by preprocessing the photodetector output. With a gain control circuit 
or limiter at the input to the PLL, the synchronization performance is shown to 
improve with increasing signal power and decreasing loop bandwidth. Equation (10) 
shows that the requirement for preprocessing can be removed if the transmitted 
pulse shape does not have a vanishing Fourier transform at the slot frequency. 
In which case the power spectrum of the detector output would contain a component 
at the slot frequency that can be tracked out by PLLs. However, the preprocessing 
approach is preferred since non-square pulses usually require higher transmitter 
bandwidths and higher peak powers for the sene number of photons transmitted. 

The performance of the PLL synchronized system was studied and compared to 
that for the perfectly synchronized system. It. is seen that higher signal levels 
are necessary to compensate for the imperfect synchronization due to the PLL. 
However, the loss in signal power is less than 0.1 decibel for loop bandwidths 
less than 0.12 of the slot frequency. This relatively small loss suggests that 
phase-locked loop synchronization PPM systems can be used to achieve reliable 
communications at a small increase in signal power. 
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APPENDIX A: POWER SPECTRUM OF THE PHOTODETECTOR OUTPUT 


The output of the photodetector can be modeled as a filtered Poisson process, 
i(t) -l h(t - T.) (A.l) 

(V J 

where h(t) is the impulse response of the detector filter. In general, the output 
of the photodetector defined above is not a stationary process. In fact, it is 
known that for PPM signaling, the output of the photodetector is cyclostationary 
with period T^, where T^ is the word period. Consequently, the power spectrum 
of the photodetector output cannot be defined as the Fourier transform of its 
autocorrelation function. Nevertheless, the power spectrum remains a useful 
concept when discussing the frequency response of the system. We shall therefore 
define the time-averaged power spectrum for a nonstationary process as [l,i0] 


suo EI|I T U) 2 |] 

where 

T 

I_U) - / . 

1 -T 


(A. 2) 


Substituting the definition of I t (<d) into S(j), we can rewrite (A. 2) as 
S(u>) - xt-if El / 1 ( t )i(T)e" ju,(t " T) dtdr] 


it! if E X [/ E[i(t)i(T)|x]e" ju(t " x) dtdTj 


(A.3) 


n taking the joint expectation of the photodetector output, we first take 
the expectation conditioned on the detector photocount rate l(t), and then remove 
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the condition by taking the expectation with respect to X. The conditional 
expectation of the filter Poisson process can be derived by differentiating the 
joint characteristic function $(u>j,'») 2 )» g* ven by 

♦Uj.» 2 ) - exp[J X(E)[e~ j ^ lh(t " 5) + “ 2h(T " ?) -l)dE] . 


The result is 

E[i(t)i(T)|xl 


3 2 * 


3^2 


0 


- / X(5)h(t - C)h(x - £)d£ + /X(5)\(n)h(t - £)h( x - n)d£dn . 

(A.4) 


By substituting Ec •• tion (A.4) into (A. 3) and carrying our :he Fourier trans- 
form, the • ower spectium of the photodetector output can be written as 


S(u>) - |H( w) | 2 E[ A (0) + | A(oj) j 2 ] 

T+« 2T 


(A. 5) 


where A^(a) ■ X(t)e ^ wt dt is the Fourier transform of the received photocount 
rate. For the PPM signaling scheme, A^(uj) is given by 


/^(u) - 2T<5i£i£) x b + X s P'u) l e J “ (kT « +Ck - Ts) 


(A.6) 


where X s> and C^’s are defined in Equation (3), and P( m) is the Fourier trans- 
form of the pulse shape p(t). The expectation in Equation (A. 5) is now taken 
with respect to the information sequence C . Using the sum rule I = I + I and 

i.j i-J i*j 

the fact that C^'s are independent for nonoverlapping codewords, the expectation 
in Equation (A.6) can be evaluated. The resulting power spectrum is 
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S(a) - |HU)| 2 . {(A +S 1 ) + 2 »S(u)(A 2 + + y~ |p(u)| 2 (l - | R( u) | 2 ) 

w 


+ I 1 A s 2 • |pU>| 2 l S(u - 2irk/T R ) } 
w 


(A. 7) 


where 


M-l 


R(u) - E[e“ jCkuTs ] = ^ [ e j£a)Ts 

M 4=0 


(A.8) 


and we have used Che fact that 


£ 2T C*> 


2ir6(ii)) 


(A. 9) 
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APPENDIX B: POWER SPECTRUM OF THE PREPROCESSED SIGNAL 

The output of the signal preprocessing circuit is the square of the filtered 
Poisson process 

i 2 (t) - (I Mt - O) 2 0>.l> 

where h(t) is the combined impulse response of the preprocessing filter and the 
photodetector response function. 

By substituting Equation (B.l) into (A. 3) and taking the expectation of the 
filtered Poisson process, the power spectrum of the preprocesed signal can be 
written as 

S(oj) - Um E^[ / e T ^{ r X(x)h 2 (t - x)h 2 (x - x)dx 

T>» 

+ (x(t)*h 2 (t)) (x( x)*h 2 ( t) ) + / X(x)X(y)h(t - x)h(t - y)h(x - x)h(x - y)dxdy 
+ 2(x(t)*b t) ) / X(x)h( t - x)h 2 ( T - x)dx + 2(x( x)*h( x) ) / X(y)h( - - y)h 2 (t - y)dy 
+ 4(x(t)*h(t) ) (X( x)*h( x) ) / X(x)h(t - x)h(x - x)dx + (X(t)*h 2 (t)) (X(x)*h(x)) 2 
+ (X(x)*h 2 ( x) ) (X(t)*h(t)) 2 + (X(t)*h(t)) 2 (X( x)*h( x) ) 2 }dtdx] . (B.2) 

Here the expectation is taken with respect to the photocount rate X(t). 

Inspection of Equation (B.2) shows the complexity of the power spectrum. Notice 

4 3 

that only terms involving X and X are significant, since all terms will be 
negligible in the limit of a large photocount. Equation (B.2) can therefore be 
simplified by dropping terms corresponding to lower order X's. Carrying out the 
Fourier transform on the remaining terms, the power spectrum can be written as 
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su) ■ tI-It EUa t < “ > V“ ) ’ (V“ )H( “ ) * A t ( u)h( “ ) )* 

+ 2A^.(m)H 2 (u)) (A^(a))H(u))*A^,(aj)H(co) ) 

+ 4 / ^(u' ) A^(u") A^(u' - (i)")H(w’)H(u - »' )H?a")H?u> - u^duMc, f 
+ |A T (o»)H(a»)*A T (o»)H( U )| 2 ] 

where 

H 2 (a)) = / h 2 (t)e _jut dt . (B.3) 

The expectation is now taken with respect to A^Cm), the Fourier transform of the 
photocount rate. By substituting Equation (A. 6) into Equation (B.3) and taking 
the expectation of the variables C^'s, the power spectrum of the preprocessed 
signal can be obtained. This power spectrum is very complicated and further 
simplifications are needed. First, it is assumed that the preprocessing filter 
blocks the dc component of the signal, which is due largely to the uniform back- 
ground count rate and contains no timing information. Next, note that after 
taking the expectation with respect to C^, the expression of the power spectrum 
contains integrals of thr following form 

I - / | R( tu * ) | 2 • G( u> * ) du> * (B.4) 

where R(u) is given by (A. 7), and G(o>') is some function of H and P, the Fourier 
transforms of h(t) and p(t). A plot of the function R(oj) shows that the magni- 
tude of R is appreciable only in the vicinity of a) = 2irk/T g . Therefore, if the 
function G(tu) is slowly varying over the region where R is appreciable, we m?.y 
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2 

approximate the function |R(o>) | under the Integral sign by a train of delta 
functions, 


| R( u>) | 2 * I 1 l 5(u> - 2*k/T s ) 


(B.5) 


2ir/T s 

o 2 • 2 

where T /2x - /|R(oj)| du is the area under each peak of |R(u)| • The width of 

“ 0 2 

each peak of |R(w)| is the order of o> a /M, where M is the order of the PPM and 
ui - 2x/T is the slot frequency. For most cases of Interest, H(tu) and P(u>) are 

S 8 

slowly varying with respect to R such that the approximations holds. Finally, 

3 4 

since X g » X^ for most cases of Interest, only terms involving X g and X g in the 
power spectrum will be retained. 


With these assumptions, the power spectrum of the photodetector output can 
be approximated. The result, after lengthy calculation, is given by 


S(u) - |F(u)H(u)*P(u)H(a») | 2 • l 6(u» - 2xk/T ) 

‘If m 4 I. S 


X 3 

+ — ~T / ?(“>' ) P^ uj” ) P( u> * - uDHCu'JHCw - - aj” )daj , da»" • 

(2w) Z w 


(B.6) 



46 


REFERENCES 

[1] R. Gagllardi and S. Karp, "M-ary Poisson Detection and Optical Communica- 
tions, " IEEE Trans. Commun. Technol., Vol. CT-17, No. 2, Apr. 1969, pp. 208. 

[2] S. Karp and R. Gagllardi, "The Design of PPM Optical Communication Systems,’* 
IEEE Trans. Commun. Technol., Vol. CT-17, Dec. 1969, pp. 670. 

[3] R. Gagllardi, "The Effect of Timing Errors in Optical Digital Systems," 
IEEE Trans. Commun. Technol., Vol. CT-20, No. 2, Apr. 1972. 

[4] W. C. Lindsey, Synchronization Systems in Communcation and Control . 

Englewood Cliff, NJ: Prentice-Hall, 1973. 

[5] A. J. Viterbi, Principles of Coherent Communcations . New York: McGraw- 

Hill, 1966. 

[6] R. Gagllardi and M. Haney, "Optical Synchronization Phase-Locking with 
Shot Noise Process," USCEE Report 396, Aug. 1970. 

[7] R. H. Forrester and D. L. Snyder, "Phase Tracking Performance of Direct- 
Detection Optical Receivers, 1 IEEE Trans. Commun., Vol. COM-21, pp. 1037- 
1039, Sept. 1973. 

[8] U. Mengal.l, "Tracking Properties of Phase Locked Loops in Optical Communi- 
cation Systems," IEEE Trans. Commun., Vol. COM-28, No. 9, Sept. 1978. 

[9] C. N. Georghlades, "Joint Baud and Frame Synchronization in Direct Detec- 
tion Optical Communications,” IEEE Trans. Commun., Vol. COM-33, pp. 357-361, 
Apr. 1985. 


47 


[10] C. N. Georghiades and D. L. Snyder, "Locating data frames in direct detec- 
tion optical communication systems," IEEE Tr. ns. Commun. , Vol. COM-32, 
pp. 118-123, Feb. 1934. 

[11] S. D. Personick, "Receiver Design for Digital Fiber Optic Communication 
Systems, I," Bell Syst. Tech. J., Vol. 52, July 1973, pp. 843. 

[12] D. L. Snyder, Random Point Processes . New York: Wiley, 1975. 

[13] W. A. Gardner adn L. E. Franks, "Characterization of Cyclo- Stationary 
Random Signal Processes," IEEE Trans. Inform. Theory, IT-21, January 1975. 

114] W. K. Pratt, Laser Communications Systems . New York: John Wiley and Sons, 

1969. 



CUMULATIVE LIST OF RADIO RESEARCH LABORATORY 


AND ELECTRO-OPTIC SYSTEMS LABORATORY REPORTS 
PREPARED UNDER NASA GRANT NSG-5049 

RRL Rep. No. 469 - Gard n **r, C. S. and N. N. Rao (December 1975), 
The Effects of Random Path Fluctuations on the Accuracy of 
Laser Ranging Systems. 

RRL Rep. No. 471 - Zanter, D. L., C. S. Gardner and N. N. Rao 
(January 1976), The Effects of Atmospheric Refraction on 
The Accuracy of Laser Ranging Systems. 

RRL Rep. No. 477 - Gardner, C. S. and J. R. Rowlett (November 

1976), Atmospheric Refraction Errors in Laser Ranging Data. 

RRL Rep. No. '78 - Hendrickson, B. E. and C. S. Gardner 

(December 1976), Correction of Laser Ranging Data for 
the Effects of Horizontal Refractivity Gradients. 

RRL Rep. No. 481 - Gardner, C. S. (February 1977), Statistics 
of the Residual Refraction Errors in Laser Ranging Data. 

RRL Rep. No. 486 - Gardner, C. S. (July 1977), Comparison 

Between the Refraction Error Covariance Model and Ray 
Tracing. 

RRL Rep. No. 488 - Gardnei , C. S. (September 1977), Speckle 
Noise in Satellite based Lidar Systems. 

RRL Rep. No. ' ) 5 - Gardner, C. S. and G. S. Mecherle (April 

1978), Speckle Noise in Direct-Detection Lidar Systems. 



49 

RRL Rep. No. 496 - Gardner, C. S. and A. M. Saleh (October 
1978), Speckle Noise in Differential Absorption Lida’- 
Systems . 

RRL Rep. No. 499 - Gardner, C. S. (January 1979), A Technique 
for Remotely Measuring Surface Pressure from a Satellite 
Using a Multicolor Laser Ranging System. 

RRL Rep. No. 502 - Palluch, E., J. D. Shelton and C. S. Gardner 
(May 1979), Operating Manual for the RRL 8 Channel Data 
Logger. 

RRL Rep. No. 505 - Gardner, C. S. and R. Axford, Jr. (March 

1980), Regression Models for Multicolor Satellite Laser 
Ranging. 

RRL Rep. No. 510 - Gardner, C. S. (April 1981), Analysis of 
Target Signatures for Laser Altimeters. 

RRL Rep. No. 511 - Gardner, C. S. (June 1981), Atmospheric 
Refraction Effects in Air Borne Laser Ranging. 

RRL Rep. No. 514 - Tsai, B. and C. S. Gardner (December 1981), 

Remote Sensing of Sea State by Laser Altimeters 

RRL Rep. No. 518 - Gardner, C. S. (August 1982), Optical 
Communications . 

RRL Rep. No. 519 - Ira, K. E. and C. S. Cardner (September 1982), 
Atmospheric Refraction Effects on Baseline Error in 
Satellite Laser Ranging Systems. 

RRL Rep. No. 526 - Ira, K. E., B. M. Tsai and C. S. Gardner 

(September 1983), Analysis of Short Pulse Laser Altimetry 
Data Obtained Over Horizontal Path. 


i 



RRJ, Rep. No. 527 - Tsai, B. M. and C. S. Gardner (March 1984), 
Theoretical and Experimental Analysis of Laser Altimeters 
for Barometric Measurements Over the 0ce3n. 

EOSL Rep. No. 84-001 - Lafaw, D. A. and C. S. Gardner (August 

1984) , Timing Performance of Phase-Locked Loops in Optical 
Pu-.se Position Modulation Communication Systems. 

EOSL Rep. No. 85-002 - Im, K. E. anc C. S. Gardner (April 1985), 
Estimation of the Differential Pulse Propagation Times in 
Two-Color Laser Ranging Systems. 

EOSL Rep. No. 85-003 - Chen, C. C. and C. S. Gardner (May 1985), 
Phase-Locked Loop Synchronization for Direct 1 u 

Optical PPM Communication Systems. 



